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ABSTRACT 

We compare QCD predictions with experimental data on inclusive multiparticle 
observables whereby we discuss various realizations of the idea of parton hadron 
duality which is related to a soft confinement mechanism. Special emphasis is given 
on effects from the running dg, the soft giuon interference, scaling properties and 
the role of the hmiting scale Qq in cascade evolution. 

1. Introduction 

There is a hierarcliy of Perturbative QCD predictions depending on the sensi- 
tivity to the hadronization process which is not yet understood at a satisfactory level 
within the theory. 

First there is the calculation of quantities in Perturbation theory in fixed order 
of as- These studies strongly support Perturbative QCD as the theory of strong 
interactions: The measurement of the total cross sections in e^e~ or in deep inelastic 
lepton nucleon scattering, then the analysis of jet cross sections in hard collisions. 
In the latter application a (not entirely trivial) assumption has to be made on the 
equality of cross sections for production of hadron jets and of parton jets at the same 
resolution. The spectacular successes of such predictions have established the basic 
properties of the coupling and the vertices in PQCD at large in terms of one scale 
parameter A;^. 

Secondly, by applying the leading log approximation for the multiparton pro- 
duction, one can push PQCD further and resolve the more detailed intrinsic struc- 
ture of jets.^~^ An additional cutoff parameter Qq is introduced which regularizes the 
collinear and infrared singularities of the gluon bremsstrahlung processes. Predictions 
can be obtained on various observables of a multiparticle final state. There are two 
different approaches towards comparison of such calculations with experiment. 

1. Parton- Hadron- Duality (PHD) in various forms. 
The general idea is that a parton jet resembles in some aspects a hadron jet. An 
initial argument for such a behaviour was the proof of "preconfinement" , i.e. the 
parton cascade already in the perturbative region prepares color singlet clusters with 
finite mass independent of energy. Q However, these clusters are too heavy for realistic 
phenomenological applications. So one has to assume some kind of soft hadroniza- 
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tion rnenhanism which allows to compare properties of hadronic and partonic final 
states-oH In the application one can distinguish further 

a) infrared and coUinear safe observables 

In this case the value of the observable doesn't change if a soft particle is added or if 
one particle is split into two coUinear particles. Such observables are independent of 
the cutoff Qo and have therefore a chance not to depend on the final stage of the jet 
evolution. Quantities of this type are energy flows and -correlations, global quantities 
like thrust etc. (for a review, see ref. 8). 

b) infrared sensitive observables 

Here particles are counted. Examples are multiplicities, inclusive spectra and corre- 
lations which are discussed in particular at a multiparticle conference. These observ- 
ables depend explicitly on the cutoff Qo and the QCD results cannot be compared to 
data immediately in a meaningful way. 

For such quantities one may follow two strategies: 
bl) One constructs again infrared safe quantities. This is possible if the Qo depen- 
dence factorizes so that it dops out after proper normalization. Another possibility 
is rescaling of variables which may lead to asymptotically safe quantities. 
b2) A more progressive strategy would be to interpret Qo as hadron mass (say 
Qo ~ itt-tt) and to compare the observables for a parton jet evolved down to hadronic 
scales directly to the experimental data . This procedure has 
momentum spectra ("Local Parton Hadron Duality" (LPHD) 

2. Parton cascade with hadronization model. 
The parton cascade evolves down to a cutoff Qq > rrih, specific models are then 
introduced to describe the hadronization process. Most popular are the string model 
by the Lund group! {Qq ~ 1 GeV) and the cluster model by Marchesini and WebberS 
(Qo ~ 0.3 GeV). 

In this review we discuss the infrared sensitive quantities in multiparticle pro- 
duction. Also for such observables analytical results are important and indispensable 
for a deeper understanding of the QCD predictions; only then can we discuss power 
laws, the consequences of running and scaling properties. 

The QCD calculations can be carried out analytically in the simplest case in 
the double logarithmic approximation (DLA) which takes into account the leading 
contributions from the infrared and coUinear divergencies and provides the asymptotic 
result at very high energies. The nonleading corrections are typically large and can be 
given by an expansion of the type 1 + a^/a^ + bag + ■ ■ ■ Such Modified LLA (MLLA) 
with sub leading corrections take into account, for example, energy recoil effects but 
not yet angular recoils (see, for example, ref. 6). Alternatively one may derive QCD 
results with Monte Carlo methods which are typically more accurate and can serve 
as an important check of the approximations of the analytical calculations. 

The aim of the study of multiparticle phenomena in this context is to learn 
about the soft limit of QCD and color confinement, the part of the theory not well 
understood. In particular we are interested how the running of a^, the soft gluon 
interference ("angular ordering"), various scaling properties, the limiting scale Qo 
intrinsic to the parton cascade are reflected in the observable particle distributions. 
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Figure 1: Average charged particle multiplicity for different CM energies ^/s and NLO QCD fit 
for running Also shown is a fit with fixed (Xg (power law) and the effect of disregarding the 

soft gluon interference 



2. Particle Multiplicities 

2.1 Total Multiplicity 

The multiplicity of partons emitted from a primary parton of momentum P 
into a cone of half opening angle behaves asymptotically in QCD like 



/P0x27o{pe) 



For the full event one sets = 1. Here 7o = y Otts/vr is the QCD anomalous dimension 
controlling the multiplicity evolution in DLA.El For running as we have 7o(pt) = 
(3^/\n{pT/A) with = U^fN^ - |A^/)"^. The prefactor comes from the NLO of 
the anomalous dimensionEl and b = 0.4916 for Nf = 5. For fixed as the multiplicity 
behaves as in (|l]) but with constant exponent 70 and rises like a power with energy. 
For running the rise is more slowly. The Qo dependence sits in the prefactor, 
therefore the ratio n{P)/n{Po) is infrared safe. Results on n of higher accuracy and 
for the full energy range are also available.! 

A fit of type (||) is shown in Fig. 10 The distinction between fixed and 
running as is not possible with present data but the difference will become larger at 
higher energies. The slope is well reproduced by QCdI with soft gluon interference 
taken into account. Neglecting the interference would increase the slope by \/2 clearly 
inconsistent with data for any reasonable A parameter. 



tFor A — 0.145 corresponding to ^^(Af^) = 0.117; a change of A by a factor of 10 would change 
the slope by ^ 15%.) 
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Figure 2: The KNO multiplicity distribution in a; = n/n for infinite energies (thin line) and LEP 
energies (thick line), calculated for a gluon jet with fixed in NNLO of QCD.El The negative 
binomial distribution (open points) with parameter K = 7 is also shown for comparison. 

2.2 Multiplicity Distribution 

In the DLA, valid for asymptotic energies, one deriveij a scaling; property for 
the probability -P„ in the rescaled multiplicity n/n ("KN0-scaling"ll3llZl). 

nPn = f{n/n) (2) 

There are large corrections to the asymptotic behavior from momentum conservation 
alone which has been studied analyticallylil0 and numerically^ In Fig. 2 we show 
as an example a calculation for a gluon jet with fixed in next-to-next-to leading 
order .Ej The asymptotic behavior for large x is exponential f{x) ~ e~^°^ (with 
j3o ~ —2.552). The approach to this scaling limit is very slow and the preasymptotic 
distribution at LEP energies looks quite different, at large x the distribution /(x) 
drops faster, like 

f{x)^exp{-[DxY), /i = (l-7)-'>l, D^C{2^/Txy (3) 

where the anomalous dimension 7 0.41 at LEP energies in this approximation and 
C ^ 2.5527. 

Analytic results for the realistic case of quark jets with running ag have been 
derived for moments of the multiplicity distributions Rq =< n{n — 1) . . . (n — g + 1) > 



5 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 



1 1 

- leading 


1 

order 


1 

— 


QCD 






next— to- 


-leading 


order 


^li— ^ 7 







Lund 7.2^ 






• 


ALEPH 






DELPHI 




■ 


JADE 




_ X parton MC o 
1 1 


TASSO 
1 


1 



20 40 60 80 100 
(s [GeV] 



Figure 3: Two particle multiplicity moment i?2 =< n{n — 1) > / < n >^ vs. CM energy. The 
QCD results in leading order (DLA), NLO0 and from parton MC (HERWIG) in comparison with 
experimental data.0 



/ < n >^. These normalized quantities are again infrared safe. The QCD prediction 
on i?2 in comparison with experimental data are shown in Fig^ 3. The DLA result at 
infinite energies is i?2 = The next-to- leading order resultcil 
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R2 = ^(1 - xVal) 



(4) 



with X = 0.55 reduces the moment by ~ 30% but still differs from data by ~ 10%. We 
have calculated this moment also with the MC method (using the program HERWIG0 
at the parton level!) which fully takes into account energy momentum conservation. 
This result finally matches the data and confirms PHD for an infrared safe quantity. 



3. Inclusive One-Particle Distributions and the LPHD Hypothesis 

3.1 Momentum Spectra 

One defines the logarithmic variable ^ = ln(P/A;) = ln(l/x) for a particle of 
momentum k emitted from a primary parton P; < ^ < Y with Y = ln{PQ/Qo). 
The asymptotic DLA results are known for fixed a^i" 



dn 



1/2 



Y 



exp(-27o(e--)V>^) 



(5) 
(6) 



§We take the perturbative cascade without g ^ qq splitting at the end of the cascade and choose 



parameters A = 0.15 GeV, m„ 



0.32 GeV. 
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and for running as where one obtains in Gaussian approximatio 

* ^ ? exp ( 'MZDLJ] (T) 

d( ((y + A)3/2 - A3/2)3 ^ (y + A)3/2-AS/2; 

with A = ln(Qo/A), P defined above. The constant as result can be recovered from 
this by the formal limit /?, A — oo, = 70 fixed. This Gaussian shape, known 

as "hump-backed plateau" , is characteristic of the destructive soft gluon interference 
resulting in the suppression of small momenta. The maximum occurs for ^* = ^. 
The width behaves like 

cr^ ~ F for a, fixed (8) 



. 3 



a r^Y2 for as running (9) 

So the running of as influcences the shape of the distribution. As Y depends on Qq 
this distribution is not infrared safe. However one can consider a high energy limit 
using the rescaled variableO 

^ (F + A) ln(P0/A) ^ ^ 



In great analogy to the case of angular correlationg^s to be discussed below we can 
construct a quantity with a scaling limit and obtain from (|^) 

InrfnM ^^ 3(C-Co/2)^ ^ 1_3 

Inn 2(l-(l-Co)^) ^ 2^ ^ ^ 

where the last limit holds for large Y , (q = Y/ {Y + X) ^ 1. So this quantity becomes 
infrared safe for high energies. 

Results in higher orders have been obtained by the St. Petersburg groupliil (for 
a review of various approximations, see ref. 6). Again the corrections to the DLA are 
sizable. For example, the position of the maximum ^* is shifted away by an amount 
of 0(^/51) 



5. = r(i + a.Hll-a^^ + ...). a = ilA'. + 2^ (12) 

^ 4 V32A',7r 32N,7r ^' 3 3^-2 v ; 

In the application to experiment these authors proposed the hypothesis of "Local 
Parton Hadron Duality" 



dn dn 
— = const— 

d^ hadron 



(13) 

parton,Qo=mi^ 



a proportionality of the spectra of hadrons to the one of partons for a cascade evolved 
down to a hadronic mass, Qq = rrih, which is taken typically as mT^. Of course, such 
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Figure 4: Distributions in ^ = Inl/Xp of charged hadrons at different energies compared with 
the analytical MLLA formula and a distorted Gaussian for parton cascade evolved down to a mass 
of 250 MeV, taken from RefS 



a hypothesis cannot be correct for all observables (for example, the mass spectrum 
of two partons does not show a p-meson), but it is interesting to explore its validity 
for sufficiently inclusive quantities. The relation (13) is tested in Fig. 4, where on 
the r.h.s. the "limiting spectrum" in MLLA is inserted which is obtained from the 
general formula by letting Qq = nth —>■ A. For charged particle spectra an effective 
hadron mass of rrih = 253 MeV is taken. It is remarkable how well the shape of the 
distribution and its energy evolution is fitted by the theoretical prediction in terms 
of 2 parameters {mh, const) for not too small momenta [k ^ 400 MeV). A closer 
inspection of Fig. 4 also shows that the deviation from the DLA prediction (|^ both 
in the position of the peak ^* and in the energy dependence are well supported.! 

Recently it has been questionedH as to what extent this result provides evi- 
dence for the soft gluon interference of QCD. A model has been considered (JETSET 
7.3) with a parton cascade, cut off at Qo ~ 1 GeV, either including or not including 
the soft gluon coherence, followed by string hadronization. The parameters of the 
hadronization model are adjusted in both cases so as to reproduce the main features 
of the data. In Fig. 5 one can see the ^-spectra of gluons, from the coherent calcu- 
lation with Gaussian shape and the non-coherent calculation with a higher particle 
density towards larger ^, i.e. smaller momenta. Both models after hadronization yield 
the same spectra for charged particles which are also in good agreement with data on 
charged particles of Fig. 4. Therefore within this hadronization model, as stated,cil 
no evidence can be claimed from the measurements in favor of the soft gluon inter- 
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Figure 5: Distributions in ^ = In 1/Xp from the string model (JETSET 7.3) at the parton level 
("Gluons") with and without the color coherence of the QCD included, and for charged particles 
(C^) of the final hadronic state (from ref.24). 



ference effect. On the other hand, evolving the parton cascade further according to 
QCD from Qo ~ 1 GeV to Qo = ^h, only the gluon distribution for the coherent 
case would approach the data (as is asserted by Fig. 4). This example nicely demon- 
strates the predictive power of the LPHD hypothesis and the similarity of parton and 
hadron spectra at comparable scales. On the other hand hadronization models have 
considerable flexibility to bring quite different theoretical schemes at the parton level 
into agreement with the data. In particular, the LPHD relation (|13]) is a consequence 
of the string model only for a special set of model parameters. 

3.2 Hadron mass effects 

In a more speculative application of the LPHD hypothesis the distribution of 
heavier particles of mass M is given as in Eq. (0) with the cutoff Qq = M > A, 
whereas for pions Qq = m.,^ ^ A. In the DLA the peak position = so for a 
heavier particle the peak is shifted by 

1 M 

A?=--l.- (14) 

with respect to the pion independent of energy. The numerical studies of the MLLA 
equationBa conform this energy independence and = /(M/A). In Fig. 6a one can 
see that the data for 7r°, rj are indeed separated by an energy independent amount. The 
peak position ^* at LEP-energies for various hadrons are shown in Fig. 6b together 
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Figure 6: Peak position ^* of the distribution in ^ = ln(P/ k) vs. CM energy for 7r° and rp3 (left 
side); ^* vs. hadron mass M at LEP energies compared with predictions from DLA and MLLA. 



with the prediction from DLA, Eq. (|T^ and MLLA extracted from the pubhshed 
results.il The ^* value drops when going from vr to K/rj similarly to the MLLA 
prediction but then saturates: the relation Qq = M does not work well for heavier 
particles. 

4. Two Particle Correlations 

4.1 Momentum Correlations 

The normalized correlation function is defined by 



in terms of the n-particle density p^^\ It has been calculated in NLO as an expansion 
in the arguments up to second orderEl and is found to depend only on the rescaled 
observables Q = C,i/(Y + A) at high energies. Alternatively one could choose the 
rescaled observable ^/^* as ^* = F/2 in DLA. In Fig. 7 we show -R2(^i, C2) for ^1 = ,^2 
vs. ^/C,* (where ^ = ^1) at the parton and hadron level as we obtained from the 
HERWIG MC at two different primary quark energies. 

Contrary to the integrated moment F2 (see Fig. 3) PHD does not work well 
for the differential moment but it improves with increasing energy, so one may speak 
of asymptotic PHD. 

4.2 Azimuthal angle correlations 

Energy-multiplicity- multiplicity correlations in the aximuthal angle around the 
jet direction have been calculated in dlaS and next-to-leading order.EJ The NLO 
calculations for large relative angles [cp = vr) differ from the MC results at the parton 
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Figure 7: Two-particle momentum correlation function R2{^i, ^2) for ^1 = ^2 vs. ^/^*, where 
^ = ^1, and ^* is the peak position of p^^\ for hadrons and partons from the HERWIG MC together 
with the analytical result in linear approximation by Fong and Webber (A=255 MeV). Also shown 
are the OPAL data.H 



level by ~ 10% but the experimental data follow the MC results closely for not 
too small relative angles in favor of PHD. The difference between coherent and 
incoherent models is rather small, typically about 5%. 

4.3 Polar angle correlations 

The correlations in the relative polar angle -^u of twonartons within the 
forward cone of half opening angle 6 have been derived in DLA.e21e3 This is a special 
case of n-particle angular observables h^'^^{S,'d, P), like multiplicity moments, to be 
discussed below. For such quantities the leading asymptotic behavior is given by 

h'^''\6,^,P) ~ exp{2(3^ln{P^/A)u;{e,n)) (16) 

. = (17) 

Because of 5 ^ Qo/P ~ A/P the rescaled angular variable fulfils < e < 1. The 
scaling function a;(e, n) is known^l and can be expanded for small e like 

1 — 1 

Lj(e,n) = n e + ... (18) 

2 n 

and for large n like 

uj(€,n) =ny/l-£ + 0(-) (19) 

n 

For the 2-particle correlation within the cone of half opening the quantity 

fi'du,e,P) = p^'\^i2,e,P)/n\e,p) (20) 
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is considered where n is the multiphcity in the forward cone. In DLA one obtains the 
asymptotic prediction 



f(^9i2, 0, P) ~ exp [2(3^\n{PQIK){uj{e, 2) - 2) j (21) 

An interesting property of this result is the "e-scahng" . Up to a known prefactor in 
the exponent the correlation function depends on the three variables 79i2,0,-P only 
through the single rescaled angular variable e. This scaling property is checked by the 
MC calculation at the parton level in Fig. 8 which shows the e-dependent part of the 



exponent in (pi]) vs. e. For sufficiently small e ^ 0.5 this scaling property is indeed 



well satisfied for P ^ 20 GeV. Eq. (|2T|) can also be seen to reproduce the trend of 
the MC data. Some nonasymptotic corrections within DLA different for quark and 
gluon jets, are also known.a In Fig. 8 we have adjusted the overall normalization 
because this is a nonleading effect. 



For small e with (plSf ) the correlations become asymptotically power behaved 

f(^i2) = (e/^i2)~^^°(^®) (22) 

and this can be related to the selfsimilarity of the jet cascade. This power law holds 
asymptotically in the full angular region for fixed a^, so the curvature in Fig. 8 in the 
asymptotic curve refiects the running of as- 

A comparison of the parton and hadron MC reveals that in the region with 
scaling e ^ 0.5 there is also PHD (in the MC). This region is again characterised 
by its independence of Qq in confirmation of the general rule. A preliminary result 
from the DELPHI collaboration,!! also presented at this conference, supports indeed 
PHD for this observable, i.e. the close similarity of the experimental and parton MC 
data (as Qo > rrth in the MC we talk here about PHD and not LPHD). Furthermore 
the first evidence for e-scaling in the opening angle (for 30° < O < 60°) has been 
presented. 

Another interesting aspect is the sensitivity to the soft gluon interference 
which is most naturally observed in the polar angle correlations. This is demon- 
strated by an analysis with L3 dataS considering the Particle-Particle-Cprrelations 
(PPC). They are defined like the well-known energy-energy correlations!^ but with- 
out the energy weights in the double sum over particle pairs. Likewise defined is 
the asymmetry PPCA {§12) = PPC (180° - 1^12) - PPC(^9i2). Disregarding the fiuc- 
tuations of multiplicity Nch this quantity is related to our correlation function like 

PPCA(^12) ^ -iP^'K^12)-p^corr.i^l2))/n'. 

In Fig. 9 results are shown from the JETSET MC for the angular asymmetry 
PPCA which varies strongly for small if angular ordering (AO) is properly included 
as required by QCD. The data follow closely the AO-case. The sharp dip at small 
')!}i2 is related to the peak in p^'^\'di2) ~ '^/'^i2 in DLA,il but also Bose- Einstein 
correlations could be important in this region. 

The angular ordering condition implies that a soft gluon is emitted from a 
parton 1 only within an angular cone limited by the next color connected parton 
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Figure 8: Rescaled 2-particle polar angle correlation vs. scaling variable 6 for different primary 
energies P and for fixed cone opening (with respect to the sphericity axis) from the parton 
MC (HERWIG) . The full curve represents the high energy limit of the DLA; the dashed curve the 
prediction for quark jets at 45 GeV; the normalization of the curves is adjusted. The insert shows 
the energy dependence of the same quantity for fixed E; from Ref.0 



{^12 < ''^i next)- We may estimate roughly {}i next ~ ^12 ~ ©/v^ which corresponds 
to about l/v^TSrad ~ 15°. One expects the emission within this angular region 
to be enhanced, outside to be suppressed in comparison to the not ordered case. 
Interestingly, such an effect is indeed seen in Fig. 9 (note that PPCA ~ — p^^^), so 
this interference effect between gluons is actually visible between hadrons in the region 
15-50°. This consideration also makes clear that the effect of angular ordering is most 
clearly seen in the polar angle (and not in the azimuthal angle or between momenta). 
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5. Angular Correlations of General Order n 

The general n-particle cumulant correlation function has also been studied in 
DLA.0 Here we discuss only as application the integral over the correlation function 
in certain angular regions. First we discuss the sidewise ring where the polar angles 
i)i with respect to the initial parton are within the range i) — 6< + 6, second 

the sidewise cone centered at a polar angle and half opening 6. We refer to these 
configurations as dimension D = 1 and D = 2 cases according to the phase space 
volume 6^. The factorial multiplicity moments are obtained from 

f^''\^,6) = J p'^^\Q,,...Qn)dni...dQn (23) 

or normalized F*^"-* = f^'^^/n". In analogy the cumulant moments are constructed 
from the cumulant (connected) correlation functions and are related to the factorial 
momentsta (i.e. C2 = F2 — 1, etc.). 

At high energies one obtains in DLAHeHIej 

CW(^9, 6) = {d/5Y- , = D{n - 1) - 2^^{Pd){n - u{e, n))/e (24) 

In leading approximation the difference between factorial and cumulant moments 
vanish.ll Corrections in the MLLA are also givenii and amount to typically 10%. 



^^1^ formula ( p^ ) is given for the cumulant moments in rcf.L^I and for the factorial moments in 
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In the hmit of small e (sufficiently large opening angles 5) one can use the 
linear approximation for uj{e,n) ( |T^ ) and ( p4D becomes a power law with 

ifn = D{n - 1) - (n - -hoiP^) (25) 

n 

This power behaviour reflects the fractal structure of the selfsimilar cascade. In this 
limit, the behavior of moments is also independent of Qq, i.e. infrared safe. For small 
angles S ^ Qo/P the above asymptotic formula is inappropriate and strong sensitivity 
to Qo appears. 

This type of power behavior was studied intensively in the last years in the 
context of "intermittency" .0 Whereas these phenomenological studies concentrate 
mainly on the small angle region the power behavior in QCD occurs only for the fully 
developed cascade at sufficiently large angles (small e). 

As in the case of two particle correlations we can exhibit the £-scaling by 
considering the quantity 

r.(n) _ ln[(V^)^("-^)C(")] 

in the high energy limit. One obtains 

-C^") ~ 2p{l -uj{e,n)/n) (27) 
^ 2/3(1 -/r^) (28) 

where the last approximation follows for the large n limit (|19D and is independent of 
n. 

In Fig. 10 we show the quantity — C*(^) for the ring (D = 1) from (^) for 
different energies from the parton MC. There is still considerable scale breaking at 
small e but the data points approach at high energies the DLA result (pT]). We also 
show in Fig. 10 the analogous results for factorial moments (definition as in (|26|)). In 
this case the e-scaling sets is already at low energies with an e-dependence as in (|27|). 
Note that the absolute scale as well as the difference between factorial and cumulant 
moments is of nonleading order in DLA. 

An interesting aspect of these results is the universal behavior of the observ- 
ables C^"'^ for different n, D and also of the quite different observable r^e). The 
observables r{€) and C*^^-* in Figs. 8 and 10 are determined from different parts of 
phase space but nevertheless approach the same asymptotic limit. Such a similarity is 
a characteristic property of the QCD cascade. The comparison of parton and hadron 
MC shows again that PHD works well in the model for not too large e near I.eJ 

6. Clans 

Clans are definedEl as group of particles of common ancestor. The clans are 
independently produced and therefore the number of clans in an event follow a Poisson 
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Figure 10: (a) Rescaled cumulant moments for the ring (D=l) as defined in Eq.(p6[) from the 
parton MC for different jet momenta P in comparison with the asymptotic prediction Eq.(p7|), (b) 
same as (a) but for factorial moments. 



distribution. With the assumption of a logarithmic distribution of particles in a clan 
one obtains a negative binomial distribution of particles, which describes the data in 
a good approximation. 

The analysis of multiplicity distributions in symmetric rapidity windows \y\ < 
He has revealed the interesting result that the average number of clans iV(?/c, a/s) at 
fixed CM energy y/s grows linearly with for a large range of Uc before it bends if 
He approaches the full phase space Ufps- There is only a weak energy dependence of 
this phenomenon. 

One may think of the clans as the jets which evolve from the gluons radiated 
off the primary partons in the collision in a Bremsstrahlung fashion with a Poisson 
distribution. This connection between parton and hadron level has been studied in 
hadronization models. In particular, the 1/k parameter of the negative binomial 
distribution was found about equal at parton and hadron level for ^/s > 200 GeV 
whereas the multiplicity n was larger for hadrons than for partons. It is interesting 
to note that 1/k parameter determines fully the normalized factorial moments. They 
are expected to be infrared safe as it is known for the full interval (see 2.2 and Fig. 
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3). On the other hand n is clearly Qq dependent and therefore different at parton 
and hadron level. In a generalization of these findings it has been suggested that 
the differential rapidity distributions p^^\yi ■ ■ ■ Un) at the narton and hadron level are 
proportional ("Generalized Local Parton Hadron Duality" This relation certainly 
works best if the cutoff Qq is close to the hadron mass. Otherwise there will be 
different kinematic limits. In that case there may be a better duality for the rescaled 
rapidities y/vmax- 

An interesting scaling property in terms of a rescaled rapidity has recently been 
found@ for the clan multiplicity. The multiplicity ratio tt* = N{yc-, Vs)/N{yfps, \/s) 
is calculated analytically in a simplified version of the QCD parton shower. This 
quantity scales well in the rescaled rapidity y* = yc/yjps 

n*{y:,^s)^7r*iy:) (29) 

in the studied range 50 < y/s < 500 GeV and approaches in the very high energy 
limit 

v/5).,:+o(^) (30) 

where the approach to this limit is very slow. 



7. Limitations of Parton Hadron Duality 



For the inclusive observables discussed so far the PHD concept seems to work 
rather well as seen from experimental data or suggested by MC calculations. We 
have, finally, to discuss where we expect or observe the limitations. 

We clearly expect deviations in the short range correlations due to resonance 
effects which occur at the hadron but not at the parton level, i.e. for masses ^ 
1. — 1.5 GeV. This causes a violation of PHD, for example, in the correlation function 
r{e) for e — 1 or 'di2 — >■ QcJ-P- A recent example of this type has been presented 
by the OPAL collaboration.Ea The measured ratio of sub-jet multiplicities in two- 
and three-jet events deviates from the QCD predictionc^ at the parton level below 
a resolution scale of ~ 2 GeV. As this scale is a bit large it would be desirable to 
understand better the reasons for the sudden change in slope at the hadron level 
(charm production?). In any case the discrepency between data and theory is only of 
the order of 10%. I 

Another interesting limit wher PHD may fail and is expected to fail in standard 
hadronization models is the quasielastic limit in e~^e~ annihilation with a large rapid- 
ity gap. Such events at the parton level correspond essentially to e'^e^ — > gg+(few 
soft gluons). The probability for this state is given by the Sudakov form factor for no 
radiation into this angular interval. With PHD the parton final state would transform 
into a similar final state with large rapidity gap (see Fig. 11). 

In the standard hadronization models there is a color field between the sep- 
arated partons which would decay into many hadrons of small pt- Assuming an 

II I would like to thank T. Sjostrand for bringing this result to my attention and for correspondence. 
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Figure 11: Hadronization of a quasi-exclusive partonic state yields a different hadronic final state 
according to Parton Hadron Duality and the string model. 

independent emission of particles or clusters with rapidity density c the probability 
for a gap of length Ay iscZl P = exp(— cAy), so there is this additional suppression 
to obtain events with a large rapidity gap. Then there is no PHD in the standard 
hadronization models in this limit. If an effect of this type was observed it would im- 
ply color bleaching earlier in the cascade than usually thought, some kind of exclusive 
PHD.i 

An effect of this type has actually been seen in the MC_simulation of the 2-jet 
rate for small resolution scales ?/cut at parton and hadron level.Ej For ?/cut ~ 10"*^ {Kt 
scale ~1 GeV) the rate for the hadronic 2-jet events is found about hundred times 
larger than for the partonic ones. 

8. Conclusions 

It is quite remarkable that the QCD calculations on the parton cascade match 
the experimental data and this suggests indeed a rather soft hadronization mecha- 
nisme. The scheme of PHD is well-defined and economic with the only parameter 
A for infrared safe quantities, and the additional quantity Qo otherwise. We have 
considered here mainly inclusive infrared sensitive quantities. 

1. PHD and LPHD. As a general rule, if the DLA result is independent of Qq (after 
appropriate normalization) then the hadronization corrections are small (multiplicity 
moments, angular observables for small e). There are cases where the Qo dependence 
disappears only asymptotically, then also PHD does the same (momentum corre- 
lations). In some cases the nonasymptotic corrections are large (cumulant angular 
correlations C^e)). The LPHD hypothesis with Qq = nih is successful for momentum 
spectra of the light hadrons (not well for heavier hadrons), but no other application 
has been provided so far. 

2. Rescaled variables appear naturally in the description of the QCD cascade: n/n, (, 
e,y*. Especially the recently introduced angular variable e{6, {}, P) provides a new tvoe 
of scaling predictions for the QCD cascade with two redundant variables 'd, P.l§E3 
First results presented at this meeting look promising for the '(9-independence.il 
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3. Soft gluon interference is taken into account by the angular ordering prescription. 
There is evidence from the shape of the hump backed plateau (see Figs. 4,5) and the 
multiplicity rise, but most sensitive are correlations in the relative polar angle ^912 
which are enhanced and suppressed at small and large relative angles. 

4. Running a^. The distributions in various observables are markedly different for 
a fixed or running «<, calculation, as in case of momentum spectra and in particular 
angular observables which are power behaved for fixed a^. 

In this review we have considered results from e~^e~ annihilations but inter- 
esting results on color coherence phenomena are coming now also from the hadron 
collider at FermilabE^I and at this conference from HERA.E^ 

Analytical QCD calculations on multiparticle final states can reveal interesting 
scaling and universality patterns which provide new insights into the intrinsic struc- 
tures of jets. With sufficient care in the selection of infrared safe quantities there is 
a promising path in the calculation of multiparticle phenomena from basic principles 
of QCD. Also the extension towards small scales Qo ~ rrih is worth pursuing. A 
statisfactory understanding though of why PHD or LPHD works so well is not really 
known. 
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